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Abstract 

In 1992, V. Weispfenning proved the existence of Comprehensive Grobner Bases (CGB) 
and gave an algorithm to compute one. That algorithm was not very efficient and not 
canonical. Using his suggestions, A. Montes obtained in 2002 a more efficient algo- 
rithm (DISPGB) for Discussing Parametric Grobner Bases. Inspired in its philosophy, 
V. Weispfenning defined, in 2002, how to obtain a Canonical Comprehensive Grobner 
Basis (CCGB) for parametric polynomial ideals, and provided a constructive method. 

In this paper we use Weispfenning's CCGB ideas to make substantial improvements 
on Montes DISPGB algorithm. It now includes rewriting of the discussion tree using the 
Discriminant Ideal and provides a compact and effective discussion. We also describe the 
new algorithms in the DPGB library containing the improved DISPGB as well as new 
routines to check whether a given basis is a CGB or not, and to obtain a CGB. Examples 
and tests are also provided. 

Key words: discriminant ideal, comprehensive Grobner bases, parametric polynomial system. 
MSC: 68W30, 13P10, 13F10. 

1 Introduction 

Let R = k\a] be the polynomial ring in the parameters a = ai, . . . , a m over the field k, and 
S = R[x] the polynomial ring over R in the set of variables x = x\, . . . ,x n . Let >~x denote 
a monomial order wrt the variables x, >~a a monomial order wrt the parameters a and )^xa 
the product order. The problem we deal with consist of solving and discussing parametric 
polynomial systems in S. 

Since Grobner bases were introduced various approaches have been developed for this 
problem. The most relevant ones are: 



• Comprehensive Grobner Bases (CGB) [We92]. 

• Specific Linear Algebra Tools for Parametric Linear systems [Si92]. 

• Dynamic Evaluation [Du95]. 

'Work partially supported by the Ministerio de Ciencia y Tecnologfa under project BFM2003-00368, and 
by the Generalitat de Catalunya under project 2005 SGR 00692 
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• Newton Algorithm with Branch and Prune Approach [HeMcKa97]. 

• Triangular Sets [Mor97]. 

• Specialization through Hilbert Functions [GoTrZaOO]. 

• DISPGB Algorithm [Mo02]. 

• Alternative Comprehensive Grobner Bases (ACGB) [SaSu03]. 

• Canonical Comprehensive Grobner Bases (CCGB) [We03]. 

This paper describes some improvements made on DISPGB. Trying to solve some of the 
examples given in the references cited above using the improved DISPGB has been an interesting 
challenge (see section 5). 

In [We92] , Professor Volker Weispfenning proved the existence of a Comprehensive Grobner 
Basis CGB wrt >~x for any ideal I C S such that for every specialization of the parameters 
Ga : R — > K' extended to R[x] — > cr^ (CGB) is a Grobner basis of the specialized ideal 

(Ja(I). He also provided an algorithm to compute it. There are two known implementations 
of this algorithm [Pe94, Sc91]. 

In [Mo95] and [Mo98], A. Montes used classical Grobner bases theory to study the load-flow 
problem in electrical networks. V. Weispfenning recommended him to use the Comprehensive 
Grobner Basis algorithm [We92, Pe94] for this problem. The use of CGB in the load- flow 
problem provided interesting information over the parameters, but was rather complicated 
and not very efficient. Moreover, it was not canonical, i.e. it was algorithm depending. 

In Montes [Mo02] provided a more efficient algorithm (DISPGB) to Discuss Parametric 
Grobner Bases, but it was still non-canonical. DISPGB produces a set of non-faithful, canoni- 
cally reduced Grobner bases (Grobner system) in a dichotomic discussion tree whose branches 
depend on the cancellation of some polynomials in R. The ideas in DISPGB however, inspired V. 
Weispfenning in [We02, We03] to prove the existence of a Canonical Comprehensive Grobner 
Basis (CCGB) as well as to give a method to obtain one. 

The main idea for building up the canonical tree is the obtention of an ideal J C R, 
structurally associated to the ideal I C S and the order which clearly separates the 
essential specializations not included in the generic case. Let us denote J as the Weispfenning 's 
discriminant ideal of (I, >~x)- in t ne new Weispfenning's algorithm, J must be computed at 
the beginning of the discussion using a relatively time-consuming method. The discriminant 
ideal was one of the lacks of the old DISPGB and an insufficient alternative algorithm GENCASE 
was provided. 

In this paper we obtain, following Weispfenning, a discriminant ideal denoted as N, which 
can be determined from the data obtained after building the DISPGB tree using a less time 
consuming algorithm and, moreover, we prove that J C N. We conjecture that J = N. We 
have verified it in more than twenty different examples, and no counter-example has been 
found. The ideal N allows to rewrite the tree getting a strictly better discussion. 

We also prove that for a large set of parametric polynomial ideals (at least for all prime 
ideals /) the discriminant ideal is principal and in this case we have a unique discriminant 
polynomial to distinguish the generic case from the essential specializations. All the theoretical 
results commented above are detailed in section 2. 

In section 3, we describe the improvements introduced in the algorithms. We have made 
a complete revision to the old release simplifying the algorithm and highly increasing its 
speed. New routines CANSPEC and PNDRMALFDRM which perform semi-canonical specifications 
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of specializations and reductions of polynomials are given. The algorithm has been completely 
rewritten and the flow control has been simplified. Further reductions of the tree, eliminating 
similar brother terminal vertices, have been performed using algorithm COMPACTVERT. 

Following P. Gianni [Gi87], we are interested in guessing whether some basis of / is a 
comprehensive Grobner basis or not, in particular for the reduced Grobner basis of I wrt the 
product order >~xa. We give, in section 4, a simple algorithm ISCGB which uses the DISPGB 
output tree to answer that question. We also give an algorithm PREIMAGE to compute a 
faithful pre-image of the non-faithful specialized polynomials from the reduced bases. This 
allows to construct a CGB. It will be interesting to compare our CGB with Weispfenning's 
CCGB when implemented. 

Finally, in section 5, we give two illustrative examples and a table of benchmarks for DISPGB 
applied to several parametric systems from which the power of the algorithm is clearly shown. 

It is stated in the same section that the new DISPGB 1 algorithm is efficient and provides a 
compact discussion of parametric systems of polynomial equations. An incipient version of it 
was presented in [MaMo04]. 

2 Generic Case, Discriminant Ideal and Special Cases 

Let K = k(a) be the quotient field of R and IK the ideal I extended to the coefficient field K. 
Consider G = gb(IK, >~x), tne reduced Grobner basis of IK wrt As K is a field, G can 
be computed through the ordinary Buchberger algorithm. The polynomials in G have leading 
coefficient 1. With this normalization g can have denominators in R. Let d g G R be the least 
common multiple of the denominators of g. To obtain a polynomial in S corresponding to g 
it suffices to multiply g by d g . Following Weispfenning [We02, We03], for each g G G we can 
obtain a minimal lifting of g, a g g, such that a g g G / and a g G R is minimal wrt >~a- Doing this 
for all g G G we obtain G' , a minimal lifting of G which Weispfenning calls the generic Grobner 
basis of (I, >-%). Of course, d g \ a g . We will use a sub-lifting of G, G" = {d g g : g G G} C S, 
and this will be our generic case basis because it is simpler to compute and corresponds to 
our standard form of reducing polynomials, as it will be seen in section 3. 

We call singular specialization a specialization a for which the set of lpp (leading power 
products) of the reduced Grobner basis of a(I) is not equal to the set of lpp(G, >- x). 

DISPGB builds up a binary dichotomic tree T(I, ^„) branching at the vertices whenever 
a decision about the cancellation of some p G R has been taken. Each vertex v G T contains 
the pair (G V ,Y> V ). T, v = (N V ,W V ) is the semi-canonical specification of the specializations 
in v, where N v is the radical ideal of the current assumed null conditions (from which all 
factors of polynomials in W v have been dropped), and W v is the set of irreducible polynomials 
(conveniently normalized and reduced by N v ) of the current assumed non-null conditions. 
Considering W* the multiplicatively closed set generated by W v , then G v C (W*) -1 (K[x]/N v ) 
is the reduced form of the basis of a (I) for the specification of the specializations a G £„. At 
a terminal vertex, the basis G v is the reduced Grobner basis of <t(J), up to normalization, for 
all specializations a £ T, v . 

Weispfenning [We02] introduces the following ideal associated to each g G G: 



Release 2.3 of the library DPGB, actually implemented in Maple and available at the site http://www- 
ma2.upc.edu/~montes/ 
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the second formula being computable via ordinary Grobner bases techniques. Then the radical 
of their intersection J = ^/fl^ec ^9 ls use d to distinguish the generic case in the algorithm. 
We call J the Weispfenning's discriminant ideal. A specialization a is said to be essential 
(for I, >-x) if J g C ker (a) for some g G G. 

V. Weispfenning proves the following two theorems: 
Wl: J = P| {ker(cr) : a is essential }. 

W2: Let a be an inessential specialization. Then 

(i) a(G) is defined for every g G G and lpp(o-(g), ^) = lpp(g, >-x)- 

(ii) cr(G) is the reduced Grobner basis of the ideal <r(7). 

In the DISPGB tree T(7, specializations are grouped into disjoint final cases i by 

the specification Ej, and for all specializations in Ej the reduced Grobner bases have the same 
set of lpp wrt >-x- 

Let 1 < i < k number the terminal vertices. We call singular cases the final cases for 
which lpp(Gj, >~x) lpp(C, >~x)- Let A be the set of indexes of the singular cases: 

A = {l<i<k : lpp^H^lpp^Hc)}- 

We denote V(7) the variety of 7 and I(V) the ideal of the variety V. The tree, being dichotomic, 
provides a partition of (K') m into disjoint sets of specifications, and thus 

(K') m = U f V(JV0 \ |J VH = C/ s (J U g , 

i=l \ wG^i / 

where C/ s is the set of points a G (K') m corresponding to singular specifications, i.e. 

17,(7, y*) = {a G (K') m : ajs is singular } = (J I V(7V i ) \ J V(«;) J . 

Theorem 1. 7ei us caH iV(7, ^) = I(U S ) the discriminant ideal. Then 

N(I,y T ) = P| TV,. 

This theorem allows to compute ./V from the output of BUILDTREE, i.e. the first tree 
construction in DISPGB. (See section 3). 

Proof. We prove both inclusions: 

C: /(a) = for all / G N = I(U S ) and a G U s . Thus o-„(/) = for all a G f7 s . Taking 
now a such that <Ta G Ej this implies that / G iVj. As this can be done for all i G A, it 
follows that iV C f] ieA Ni. 

D: For all / G f] ie ^ Ni and all a £ U s there exists i £ A such that 05- G Ej and, of course, 
/ G JVi. Thus <7o(/) = 0, i.e. /(a) = for all a G E/ s . Thus / G I(U S ) = N. 

□ 

Before proving the next theorem we need the following 
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Lemma 2. Any singular specialization is essential. 

Proof. Let o~a be a singular specialization. If it were not essential, by Weispfenning theo- 
rem (W2), then the reduced Grobner basis of a (I) would be the generic basis G, and this 
contradicts the definition of singular specialization. Thus a-a must be essential. □ 

Theorem 3. JCJV. 

Proof. By Weispfenning's theorem (Wl), if / G J then / G ker(<7a) for all essential o~a, an d 
thus f(a) = 0. So, by lemma 2, f(a) = for all singular ctq. This implies that f(a) = for 
all i G A and era G Ej and thus / G y^W = ^i- Finally, by proposition 1,/eJV. □ 

Conjecture 4. IVe formulate two forms 

(i) (Strong conjecture). All essential specializations are singular. 

(ii) (Weak conjecture). J ^ N. 

Proposition 5. XTie strong formulation of conjecture 4 implies the weak formulation. 

Proof. If / G iV then, for all i G A, / G iVj. Thus, f(a) = for all singular specialization a-a 
and, if the strong form of the conjecture is true, then f(a) = also for all a-a essential and 
thus / G ker(cTa)- So, by Weispfenning's theorem (Wl), / G J. □ 

In any case, by definition N is discriminant, i.e. for any a G" V(N) the Grobner basis of 
o~a(I) is generic, and every singular specification is in Y(N). Thus, what we called minimal 
singular variety in [Mo02] is described by Y(N). If the strong formulation of the conjecture is 
true then every specialization a, for which N C ker(<r), is not only essential but also singular 
and thus the corresponding set of lpp of its reduced Grobner basis cannot be generic. 

We have tested our conjecture in more than twenty examples and we have not found any 
counter-example of any of the two formulations. Nevertheless the weak formulation is the 
most interesting one and a failure of the strong formulation would not necessarily invalidate 
the weak formulation. 

In most cases Weispfenning's discriminant ideal J is principal, as states the following 

Theorem 6. If I C S is a prime ideal and the generic Grobner basis G wrt is not [1], 
then the discriminant ideal J(I, >~x) is principal and is generated by the radical of the 1cm of 
all the denominators of the polynomials in G. 

Proof. Take g G G. We have J g = d g (I : d g g) f)R. If h G J g then d g \ h, as d g g has no 
common factor with d g . Thus d g g(h/d g ) G /. By hypothesis, d g g / 1 and I is prime. So, as 
h/d g G R, we have h/d g G" /. Thus, necessarily d g g G I and d g G J g . As d g \ h for all h G J g , 

it follows that J g = (d g ) is principal. As J = ^Jf] geG J g is the intersection of principal ideals, 

the proposition follows. □ 

Not only prime ideals have principal discriminant ideals as the next example shows: Take 
/ = (ax + y + z + b,x — 1 + ay + z + b, x + y + az + b) . 
Computing the Grobner basis of / wrt lex(x, y, z, a, b) one can see that 

/ = ((a + 2)z + b,y — z,x + y + az + b) PI (a — l,x + y + az + b) 
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Routines of the 
old algorithm 


Routines of the new 
algorithm 


Improvements 


Obsolete 
routines 


DISPGB 
BRANCH 


DISPGB 
BUILDTREE 
DISCRIMINANTIDEAL 
REBUILDTREE 
COMPACTVERT 


BUILDTREE replaces old 
BRANCH 

x_>x k.ri.x» vjii • 

Current DISPGB includes 
also rpbnildinp - of t.ViP 
tree (REBUILDTREE) and 
COMPACTVERT. 


GENCASE 


BRANCH 
NEWVERTEX 


BUILDTREE 


Better flow control, no 
incompatible branching. 


BRANCH 


NEWCOND 

X \ J 1 VV \J VJ X \ LJ 


CONDTOBRANCH 


More rnf>nst ensures no 

V L I V I V i 1 U , WllV VjO _LVW 

incompatible branches. 


NEWCOND 


CANSPEC 


CANSPEC 


Uses radical ideal A/Tore 
robust. 






PNORMALFDRM 


Standard polynomial re- 
duction wrt X]. 




CONDPGB 


CONDPGB 


Uses CONDTOBRANCH and 

Weispfenning's standard 
pair selection. 






DISCRIMINANTIDEAL 


Determines the discrim- 
inant ideal N. 






REBUILDTREE 


Rebuilds the tree start- 
ing the discussion with 
N. 


GENCASE 

(external) 




COMPACTVERT 


Drops brother terminal 
vertices with same lpp 
sets. 





Table 1: 



and / is not prime. The generic Grobner basis wrt lex(x, y, z) is, in this case, G = [z + b/(a + 
2), y + b/(a + 2), x + b/(a + 2)]. Thus d g = a + 2 for each g G G. For this example it is easy 
to compute J = ((a + 2)(a — 1)) which is still principal even if / is not prime and has a prime 
component with generic Grobner basis [1]. 

It would be interesting to characterize which ideals I C S have principal discriminant 
and which do not. But it is now clear that in the most interesting cases we have princi- 
pal discriminants. This gives a new insight into our concept of singular variety used in the 
algorithm [Mo02] in order to understand the parallelism and differences between the new 
Weispfenning's algorithm [We02, We03] and DISPGB, and allows us to improve the old algo- 
rithm. 

Under that perspective, we have completely revised [Mo02] and obtained a much more 
efficient and compact discussion. An intermediate version was presented in [MaMo04]. We 
shall describe now the improvements introduced in the new DPGB library and refer to [Mo02], 
where the old DPGB is described, for all unexplained details. 
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T^DISPGB(5,Hc, >m) 
Input: 

B C i?[a][x] : basis of /, 

>-x, >s '■ termorders wrt the variables x and the parameters a respectively. 
Output: 

T: table with binary tree structure, containing (G v , £„) at vertex v 
BEGIN 
T := (/), # global variable 
v := [ ] # (label of the current vertex) 
£ := ([ ],(f>) # (current specification) 

BUILDTREE(-u, B, S) # (recursive, stores the computations in T) 
N : = DISCRIMINANTIDEAL(T) 
COMPACTVERT(T) # (compacts T) 
REBUILDTREE(T, N) # (rebuilds T) 
COMPACTVERT(T) # (compacts T) 
END 

Table 2: 

3 Improved DISPGB Algorithm 

In this section we describe the improvements introduced in DISPGB algorithm. Table 1 sum- 
marizes the basic differences between old [Mo02] and the new algorithms used in it. 

First, we have improved the construction of the discussion tree T(I, in order to 

have a simpler flow control and to make it faster by avoiding unnecessary and useless time- 
consuming computations. In the old algorithm this was done by the recursive routine BRANCH 
which was the unique action of DISPGB, but now it is done by BUILDTREE. As we explain later, 
it has been strongly reformed. 

Then, DISCRIMINANTIDEAL computes the discriminant ideal N = HieA ^» which, as shown 
in section 2, can be determined from BUILDTREE output. 

After that, DISPGB calls REBUILDTREE. This algorithm builds a new tree setting the dis- 
criminant ideal N at the top vertex and the generic case at the first non-null vertex labelled as 
[1] (see figure 1 in section 5.1). The old tree is rebuilt under the first null vertex recomputing 
the specifications and eliminating incompatible branches. The result is a drastic reduction of 
branches in the new tree. In the old DPGB library, this work was partially done by the external 
algorithm GENCASE which has become useless. 

To further compact the tree, a new algorithm COMPACTVERT is used. It summarizes brother 
terminal vertices with the same set of lpp into their father vertex. COMPACTVERT is called before 
and after REBUILDTREE. DISPGB algorithm is sketched in Table 2. 

3.1 Building up the Discussion Tree: BUILDTREE. 

We have simplified the flow control from the ancient DISPGB and dropped useless operations. 
Now all the hard work of the discussion is done by the recursive algorithm BUILDTREE which 
replaces the old BRANCH routine and makes NEWVERTEX useless. The obtained discussion is 
equivalent to the one given by the old DISPGB, but now is more compact. 

It computes the discussion tree faster than the old one because now it assembles the dis- 
cussion over the coefficients of the current basis in one single algorithm, avoiding unnecessary 
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BUILDTREE^, 5, E) 
Input : 

v, the label of the current vertex, 

B C i?[a][x], the current basis, 

E = (N, W) the current reduced specification. 
Output: No output, but the data are stored in the global tree variable T. 
BEGIN 
Cf := false 

(c b , c d , G, S , Si):=C0NDT0BRANCH(5, E) 

IF Cd THEN # is true if all lc(g), g G G are decided non-null, false 
otherwise) 

(c fe , c/, G, S , Si):=CONDPGB(G, S) 
END IF 

T v := (G, S) # (Store data in the global tree variable T) 

IF cj THEN # (cj is true if the new vertex is terminal, false otherwise) 

RETURNO 
ELSE 

IF ^ THEN # (cfe is true if null and non-null conditions are both 
compatibles) 

BUILDTREE((?;,0),G, S ) 

BUILDTREE((t>, 1),G, Si) 
ELSE 

BUILDTREE(v, G, Si) # (and BUILDTREE continues in the same vertex) a 
END IF 
END IF 
END 



"In this case, if CQNDPGB has already started then the list of known S-polynomials reducing 
to can be kept. 



Table 3: 

branching and useless computations. 

Given B, a set of polynomials generating the current ideal, BUILDTREE takes the current 
basis B v at vertex v, specialized wrt the current reduced specification S^ = (N v , W v ), builds a 
binary tree T containing the discussion under vertex v, and stores all the data at the vertices 
of T. It is a recursive algorithm and substitutes the old BRANCH and NEWVERTEX. See table 3. 

Theorem 16 in [Mo02] still applies to the reformed BUILDTREE, thus we can assert the 
correctness and finiteness of the algorithm. 

The most important algorithms used by BUILDTREE are commented below. 

The algorithm CONDTOBRANCH substitutes the old NEWCOND. It is used each time that 
BUILDTREE is recursively called and also inside CONDPGB, applying it to each new not-reducing- 
to-zero S-polynomial. This prevents Buchberger algorithm from stopping and saves incom- 
patible branches. 

Each time we need to know whether a given polynomial f E R — for example the lc (leading 
coefficient) of a new S'-polynomial- is zero or not for a given specification, we will reduce it 
by S = (iV, W) using PNORMALFORM and then test whether the remainder is compatible or not 
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(<%, Cd, G, S , Si) «- CONDTOBRANCH(B, S) 
Input : 

B C i?[a][x], the current basis 
S = (A, W) a reduced specification. 
Output : 

G is B reduced wrt S, 

Si is the reduced specification for the not null branch 
So is the reduced specification for the null branch 
Cb is true whenever So exists, and false otherwise. 

Cd is true if all g G G have lc(<?) decided to not null, and false otherwise. 
BEGIN 

G : = PN0RMALF0RM(£, S) 

IF there is g G G with l g = \c(g) not yet decided to not null wrt S THEN 
Cd := false 

(t, Si) := CANSPEC(A S , W E {J{l g }) 

(t, S ) := CANSPEC((A S , l g ), W E ) 

IF t THEN c b := true ELSE c b := false ENDIF 
ELSE 

Cd ■= true 
ENDIF 
END 



Table 4: 

with taking it null and non-null for each of the specifications using CANSPEC. The whole task 
is done by CONDTOBRANCH. See table 4. 

BUILDTREE uses a Buchberger-like algorithm -CDNDPGB (Conditional Parametric Grobner 
Basis) - taking the specification into account and intending to determine a specializing Grobner 
basis. The basic improvements on CONDPGB in the new version are: the call to CONDTOBRANCH 
instead of old NEWCOND and improving Buchberger algorithm by considering Weispfenning's 
normal strategy of pair selection [BeWe93]. We do not detail these improvements. 

CANSPEC has also been modified. 

At each vertex v of the tree a pair (G v , S„) is stored, where S^ = (N v , W v ) is a specification 
of specializations. This means that for all a G S^ one has a(N v ) = and a(w) / Vu> G 
W v . From the geometric point of view, a given S = (N, W) describes the set of points 
nN)\(U weW V(w))C(KT l - 

By proposition 5 in [Mo02], one can see that S = (N,W) and S' = (s/N,W) describe 
equivalent specialization sets. And, by proposition 7, the same happens with S = (N,W), 
where has no factor laying in W and is radical, and W is the set of the irreducible factors 
of W with multiplicity one reduced modulus N. So we choose the following representative for 
the specifications describing equivalent specialization sets: 

Definition 7. We call S = (A, W) a reduced specification of specializations if it is a specifi- 
cation such that 
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(t,£) *- CANSPEC(E) 

Input: £ = (A, W) a not necessarily reduced specification. 
Output : 
t: a boolean valued variable. 

£: a reduced specification if t = true, and otherwise (in this case 
incompatible conditions have been found). 
BEGIN 
N a := N, N b := 

W a := W, W b := the irreducible factors of W without multiplicity and 
reduced wrt N a ; 
IF Uqew b 1 = ° THEN RETURN(false,</>) ENDIF 
WHILE (N a / N b AND W a / W b ) DO 

N a :=</> 

FOR p e A fe DO 

p := drop from p all irreducible factors laying in W& 

IF p = 1 THEN RETURN(false,0) ENDIF 

Add p into N a 
END FOR 
W a := W b 
N b :=VK 

W b := {irreducible factors of W a without multiplicity and reduced wrt A^,} 

IF U q ew b 1 = ° THEN RETURN(false,^>) ENDIF 
END WHILE 
E:={N a ,W a ) 
RETURN(true, S) 
END 



Table 5: 

(i) (N) is a radical ideal, and = gb((A), >~a), 

(ii) there is no factor of any polynomials in (N) laying within W, 

(hi) W is a set of distinct irreducible polynomials not laying within (A), 
(iv) W N = W. 

We must note that the set W is not uniquely determined, as there exist infinitely many 
polynomials which cannot be null for a given specification. For example, suppose that the 
current reduced specification is W = {a}, A = [a 2 — 1]. The condition a ^ is compatible 
with A but is redundant in this case. We can also add to W other polynomials like a — 2. 
Thus there is no unique reduced specification, but our choice is convenient enough. The task 
of obtaining reduced specifications and testing compatibility of the current null and non-null 
conditions is done by the reformed CANSPEC. See table 5. 

Proposition 8. Given any specification of specializations S = (A, W), if CANSPEC (S) returns 
(t, X) with t = true, then E is a reduced specification of E computed in finitely many steps. 
Otherwise it returns t = false and (A, W) are not compatible conditions. 
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/ <- PNORMALFORM(/, E) 

Input: / 6 i?[x] a polynomial, £ = (iV, W) a reduced specification, 

Output: / reduced wrt £ 

BEGIN 

~ — N 

f := the product of the factors of / not laying in W, conveniently normalized 
END 



Table 6: 

Proof. At the end of each step iV a is a radical ideal, W a is a set of irreducible polynomials with 
multiplicity one reduced wrt N a , so W a a = W a . So, A5 is still radical when the algorithm 
stops, as Nb is built by dropping from N a all those factors laying in W a . If the algorithm 
returns true, as at each completed step (JVj,, Wb) satisfies the conditions of definition 7, then 
the conditions are compatible and S is a reduced specification of specializations. Otherwise 
the conditions are not compatible. 

Let us now see that this is done in finitely many steps. The algorithm starts with Nq = N. 
At the next step it computes Ni, and then N2, etc... These satisfy Nq C Aq C A2 C ■ ■ ■ . By 
the ACC, the process stabilizes. So, only a finite number of factors can exist, thus dropping 
factors is also a finite process. □ 

The second necessary task is to reduce a given polynomial in S wrt S. This is done in a 
standard form by PNORMALFORM. To eliminate the coefficients reducing to zero for the given 
specification it suffices to compute the remainder of the division by N, because N is radical. 
And then, in order to further simplify the polynomials, all those factors lying in W are also 
dropped from N. See table 6. 

Nevertheless, the reduction using PNORMALFORM does not guarantee that all the coeffi- 
cients of the reduced polynomial do not cancel out for any specialization a E S. To test 
whether adding a new coefficient to the null conditions is compatible with £ we need to apply 
CONDTOBRANCH. 

Given f,g £ S and £ we say that their reduced forms /s and computed by PNORMALFORM 
are equivalent wrt £ when <Ja(f) and a-a(g) are proportional polynomials for every particular 
specialization cr„ G S such that o"a(lc(/s)) 7^ and <7a(lc(<?s)) 7^ 0. 

Consider for example, £ = (N = [ab— c, ac— b, b 2 — c 2 ], W = 4>), /s = ax+c 2 , gs = cx + c 2 b 
and >~a= lex(a, b, c). /s and g% are not identical, but note that they are equivalent. As can be 
seen in this example PNORMALFORM is not always able to reduce them to the same polynomial. 
Nevertheless, we have the following 

Proposition 9. Given two polynomials f,g G S then /s ~ gs wrt E iff 

(i) lpp(/s, >~x) = lpp(#s, >~x) and 

(ii) PNORMALFORM applied to 1c(#e)/e - lc(/s)5s returns 0. 

Proof. Obviously if one of both hypothesis fail, the reduced expressions are not equivalent 
wrt S. 

On the other side, suppose that (i) and (ii) hold. Then, using order )^xa we have 

N N 

1c(#e)/e = 1c(/e)5e by hypothesis (ii). Thus, lc(# s )(a) fc(x, a) = lc(/ s )(a) g^(x, a), 
for all specializations in S. In particular it also holds for those specializations which do not 
cancel the leading coefficients of /s and g%. And so, it follows that /s and <7e are equivalent 
wrt S. □ 
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Thus, PNORMALFORM does not obtain a canonical reduction of / wrt S, but it can canonically 
recognize two equivalent reduced expressions. 

3.2 Reduction of Brother Final Cases with the Same lpp 

In many practical computations and after applying these algorithms to a number of cases, we 
have observed that some discussion trees have pairs of terminal vertices hung from the same 
father vertex with the same lpp set of their bases. As we are only interested in those bases 
having different lpp sets, then each of these brother pairs, {vq,vi}, can be merged in one 
single terminal vertex compacting them into their father v and eliminating the distinction of 
the latter condition taken in v. 

Regarding this construction, we can define a partial order relation between two trees if, in 
this way, one can be transformed into the other. 

Definition 10. Let S and T be two binary trees. We will say that S > T if 

(i) T is a subtree of S with same root and same intermediate vertices, and 

(ii) for each terminal vertex v 6 T there is in S either the same vertex v € S 
such that (G Vt , T, Vt ) = (G Vs , or a subtree S C S pending from ver- 
tex v £ S with all its terminal vertices u £ S with lpp(G u _) = lpp(G„ T ). 

So now, given a discussion binary tree T, we may find the minimal tree T within the set 
of all trees which can be compared with T regarding this relation. This is done by a recursive 
algorithm called COMPACTVERT. 

Let us just note that the minimal tree will not have any brother terminal vertices with 
the same lpp sets of their bases. 

3.3 Rewriting the Tree with the Discriminant Ideal 

The tree T built by BUILDTREE can be rebuilt using the discriminant ideal N (see section 2). 
By theorem W2, if we are given a-a such that there exists some 5 G N for which <7a(<5) / 0, 
then (Ta{I) corresponds to the generic case. Thus, placing N into the top vertex labelled as 
[ ] in the new tree T', for its non-null son vertex we will have Xj'^ = (G[i], where is 
the generic basis and £[u is a union of specifications from T corresponding to 

S[i] = {a : 35 £ N such that a{6) / 0}. 

No other intermediate vertices hang from this side of the top vertex. If the strong formulation 
of conjecture 4 holds, then no generic cases will hang from the first null vertex. 

The subtree under the top vertex hanging from the first null son, for which the choice is 
a(N) = 0, will be slightly modified from the original T. The terminal vertices corresponding 
to singular cases hanging from it will not be modified as, by construction, for all of them the 
condition is verified by the corresponding specifications. Thus we can rebuild the tree using 
the recursive algorithm REBUILDTREE which goes through the old tree T and rewrites the new 
one X". At each vertex v it tests whether the condition iV is already included in N v . If it is 
the case, then it copies the whole subtree under it. Otherwise it adds N to the null ideal N v 
and calls CANSPEC to check whether the new condition is compatible or not. If the condition 
is compatible then the basis will be reduced using PNORMALFORM and the algorithm continues. 
If it is not, then the recursion stops. This algorithm produces a better new tree with possibly 
less terminal cases (only generic type cases can be dropped). This reconstruction of the tree 
is very little time-consuming. 
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3.4 New Generalized Gaussian Elimination GGE 

We add here a short description of the improvements on the generalized Gaussian elimination 
algorithm GGE. 

We realized, by analyzing the procedure of the old GGE [Mo02], that there were some 
special cases for which we could guess the result of the divisions at each step and thus could 
be skipped. These improvements halve the computing time. 

Even though it is more efficient and faster, GGE has become not so useful now because the 
new improvements in DISPGB, detailed above, make, in general, DISPGB work faster without 
using GGE. So now, the use of GGE within the execution of DISPGB is just optional (not used by 
default). However, it can be very useful for other applications, like in the tensegrity problem 
shown in section 5, to eliminate some variables and simplify a given basis. 

4 Comprehensive Grobner Basis 

In [We02, We03] the main goal is to obtain a Comprehensive Grobner Basis. With this aim, 
we have built an algorithm, called ISCGB, to test whether a given basis G is a comprehensive 
Grobner basis for / or not. It uses PNORMALFORM algorithm to specialize G for every terminal 
case in the discussion tree. Then it checks if lpp(<r(G)) includes the set of lpp of the reduced 
Grobner basis wrt £ for every terminal case. If this is true for every final case then ISCGB 
returns true otherwise returning false. 

The algorithm also informs for which cases a given basis is not a CGB. Thus we can 
compute pre-images of the polynomials for which B does not specialize to a Grobner basis 
and add them to the given basis in order to obtain a Comprehensive Grobner Basis. 

Consider a terminal case (G v , £„) and g G G v . To simplify notations we do not consider the 
subindex v. Let H g = {/i, . . . , f r } be a basis of the ideal I g = If)(g, N) whose polynomials 
are of the form qg + n, with q £ S and n £ (N). I g contains all the polynomials in / 

which can specialize to g (for those with a(q) a non-null element of R wrt X). Set f- = fi N . 
Obviously, H' g = {/{,..., f^} is a basis of <r(I g ). Using Grobner bases techniques we can 
express g £ cr{I g ) in the form g = J2i a ifl where the c^'s are reduced wrt TV, as we are in 
Ig/N. Then h = Y2i a ifi specializes to g and is a pre-image of g in /. This is used to build 
algorithm PREIMAGE which computes a pre-image of g. 

Combining ISCGB and PREIMAGE we compute a CGB using the algorithm sketched in table 
7. Let B = gb(i~,H™), which is a tentative CGB [FoGiTrOl, Ka97], and F = {((?;,£;) : 
1 < i < k} the set of final cases of the discussion tree built up by DISPGB. ISCGB informs 
about the polynomials in F which do not have a pre-image in the current tentative CGB. CGB 
algorithm adds pre-images of them until a CGB is obtained. Nevertheless, this construction 
is not canonical and is much more time-consuming than building up the tree, because it uses 
the product order >~xa instead of working separately wrt >~x and >~a- 

5 Examples 

We have selected two significative detailed examples. The first one is the classical robot 
arm, which has a very nice geometrical interpretation, and the second one is the study of a 
tensegrity problem described by a linear system with the trivial null solution in the generic 
case which has a non principal discriminant ideal. After that, we outline a table containing 
some relevant information for several other examples. 
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B <- CGB(B,F) 
Input: 

5 = gb(J, y^g) 

F = {(Gi, Sj) : 1 < i < fc} obtained from DISPGB 
Output: B a CGB of/ 
BEGIN 
B = B 

F = SELECT cases from F for which 13006(5, ^) is not a CGB. 
WHILE F is non empty DO 

TAKE the first case (Gi, Si) G F 

5 = 5 (J {PREIMAGE(5,Si,B) : t/GGi} 

F = SELECT cases from F for which ISCGB(5, y^) is not a CGB. 
END DO 
END 



Table 7: 



5.1 Simple Robot 

The following system represents a simple robot arm (compare with [Mo02]): 

B = [s\ + c\ - 1, s% + c| - 1, / (si s 2 - ci c 2 ) - ci + r, 
Z (si c 2 + ci s 2 ) + si — z] 

Using the orders lex(si, c\, s 2 , c 2 ) and lex(r, z, Z), respectively for variables and parameters, 
DISPGB produces the following outputs: The discriminant ideal is principal: N = J = [I (z 2 + 
r 2 )]. The set of final cases expressed in the form Tj = (Gj, (iVj, Wj)) is: 

T {1] = ([2lc 2 + l 2 + l-z 2 -r 2 ,4l 2 S 2 + (l 2 -l) 2 

-2{l 2 + l)(r 2 + z 2 ) + (z 2 +r 2 ) 2 , 

2 (r 2 + z 2 ) a - 2 z I s 2 - r (r 2 + z 2 - I 2 + 1), 

2 (r 2 + z 2 ) Sl + 2 Z r s 2 + z (/ 2 - r 2 - z 2 )], ([],{/ (r 2 + z 2 )})). 
^[0,1,1,1] = ([2Zc 2 + Z 2 + l,4(Z 2 -l)r Cl + 2z/ S2 -(/ 2 -l)r, 

(Z 2 - l) 2 - 4z 2 , 4 (Z 2 - 1) z si + (Z 2 - l) 2 + 4z 2 ], 

([z 2 +r 2 ],{z,Z + l,r,M-l})), 
^[0,1,1,0] = ([1]. {[z,r],{l + 1,1,1-1})), 
%i,o,i] = ([1], ([l 2 -l,r 2 + z%{z,l})), 
^[0,1,0,0] = ([/c 2 + l, S2 ,s 2 + c 2 -l], ([Z 2 -l,z,r], {/})), 
r [0 ,o,i] = (W. (W,{r 2 + z 2 -l})), 

^[0,0,0] = ([si + c 2 -l, Cl -r, Sl -z], ([Z,r 2 + z 2 -l],{ })), 



The generic case T^j gives the usual formula for the robot. It is characterized by the 
discriminant ideal N. The singular cases have simple geometrical interpretation and give 
information about the degenerated cases. 

A graphic plot of the tree is also provided in the library. There, the deciding conditions 
can be visualized at the intermediate vertices and the lpp sets of the reduced Grobner bases 
are shown at the terminal vertices (see figure 1). 

Now we apply ISCGB to GB = gb(B,\ex(si,ci,s 2 ,c 2 ,r,z,l) wrt the output tree. The 
result is false, and the list of specializations for all the final cases is provided: 
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Figure 1: DISPGB's graphic output for the robot arm. 



[[1], {si, s 2 ci, s 2 Si, a, c 2 si, c 2 , s\ , si}, {s 1 ,c 1 ,c 2 , si}, true]] 
[[0, 1, 1, 1], {si,s 2 , s 2 ci,s 2 si,c 2 si,c 2 , sj, s^}, {si, s 2 , ci,c 2 }, false], 
[[0, 1, 1, 0], {1, si,s 2 , s 2 s 1 ,c 2 s 1 ,c 2 , s\ , si}, {1}, true], 
[[0, 1, 0, 1], {si, s 2 , s 2 ci,s 2 si,c 2 s 1 ,c 2 , s\, si}, {1}, false], 
[[0, 1, 0, 0], {s 2 ,s 2 c 1 ,s 2 s 1 ,c 2 s 1 ,c 2 ,sl, si}, {s 2 , c 2 , s\ }, true], 
[[0, 0, 1], {1, si,s 2 ci,s 2 si,c 1 ,c 2 si,s\, s\}, {1}, true], 
[[0, 0, 0},{s 1 ,s 2 c 1 ,s 2 s 1 ,c 1 ,c 2 s 1 ,sl, si}, {si,ci, si}, true], 
There are only two cases for which GB is not a CGB. Even so, the algorithm CGB only- 
needs to add one single polynomial to obtain a CGB. 

CGB = [2lc 2 + I 2 + 1 - z 2 - r 2 , c 2 2 + s 2 2 - 1, 2(z 2 + r 2 )ci - 2zls 2 
+r(l 2 — 1 — z 2 — r 2 ), Azs 2 c\ — Azrs 2 + Arc 2 c\ + Alrci 
+2{z 2 - r 2 - l)c 2 - l(z 2 + r 2 -I 2 + 3), 2rc x s 2 - 2zc x c 2 - 2zlc x 
+ (-r 2 + z 2 - 1 + l 2 )s 2 + 2zrc 2 , 2(l 2 - l)si - Alc x s 2 + 2ls 2 r 
-z(r 2 + z 2 -I 2 - 3), 2s x z + 2c x r - r 2 - z 2 + I 2 - 1, 
rsi - zc\ + /s 2 , S1C2 + ^si - cis 2 + rs 2 - zc 2 , sis 2 + c\c 2 
+/ci - zs 2 - rc 2 , c\ + s\ - 1, 4(r 2 + z 2 )c\ - 4r(l + z 2 + r 2 - / 2 )ci 
+ (r 2 + z 2 -/ 2 + l) 2 -4z 2 ]. 



5.2 Tensegrity Problem 

We study here a problem formulated by M. de Guzman and D. Orden in [GuOr04]. 

Given the five points Pi (0,0,0), P 2 (l,l,l), P 3 (0, 1,0), P 4 (1,0,0), P 5 (0,0, 1) we want to 
determine a sixth one P&{x, y, z) for which the framework with vertices {Pi, . . . , P§} and edges 
({Pi,--,p<5}^ y {PiP 6 , P 2 P 4 , P3P5} stays in general position and admits a non-null self-stress. 
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Figure 2: DISPGB graphic output for the tensegrity problem. 
The system describing this problem is the following: 

B = [Wn + W U , W12 + W 13 ,Wi2 + W 15 ,W 12 + ^23 + ^25 ~ ^26^ + ™26, 
U>12 + W 2 5 ~ W 2 %y + U>26, W\2 + ™23 ~ W2%Z + W 2 <5, W 2 3 + W 34 + XW 36 , 
Wl3 + W 3 4 - W36U + W 3 Q, ^23 + ZW 3 Q, W 14 + If 34 + W 45 - W 4 qX + W 4 Q, 

W34 + yw m , W 45 + ZW 56 ,W 15 + w 45 - zw 56 + w 56 , 

-U>26 + W26X + XW 36 - W 4 Q + W 4G X + W 5 qX, 

-W26 + W26V ~ w 3 6 + w 36 y + yw 46 + w 56 y, 

-W2G + W26Z + ZW3Q + U?46-2 ~ ^56 + ZW<r>§\ 

Set ^xa= lex(wi2,wi3,wu,wi 5 ,W23,W25,W3 4 ,w 45 ,W26,W3Q,w 4 Q,w 56 ,x,y,z). In order to sim- 
plify the system we compute GGE(B, >-xa) (Generalized Gaussian Elimination). The GGE basis 
can be expressed as B' = B[ U B' 2 , with B' 2 being the elimination ideal wrt the variables 
1^26 = W2-, W3§ = W3, w 4 q = w 4 , w^q = W5, and B[ expressing the remaining variables linearly 
in terms of W2, W3, w 4 , W5: 

B[ = [w 45 + zw 5 , w 34 + yw 4 , W25 + w 5 y, W23 + zw 3 , w 15 - 2zw 5 + w 5 , 

w 44 - 2zw 5 + W 5 , W13 - 2zw 5 + U> 5 , W12 + 2zw 5 - w 5 ] 
B' 2 = [—ZW5 + w§x — w$y, —ZW5 + w 4 z, w 4 x + yw 4 — w 4 — zw§ + W5, 

W3V -w 3 + yw 4 - 2zw 5 + w 5 , XW3 - yw 4 - zw 3 , 

w 2 z -W2 + zw 3 + 2zw 5 - u?5, w 2 y - W2 + w 5 y + 2zw 5 - w 5 , 

w 2 x -W2 + ZW3 + w 5 y + 2zw 5 - W5], 

Then, using the orders >~x= lex(u; 2 , 103, w 4 , W5) and >~a= lex(x, y, z), for variables and param- 
eters respectively, DISPGB produces the following output: 

T[i] = ([^5,^4,^3,^2], ([ ], 

{[y 2 z - yz, zx - z 2 , x 2 - y 2 - z 2 - x + y + z}})) 
r [o,i,i,i] = ([w 5 ,w 4 ,W2Z-W2 + zw 3 ],([y- l,x-z],{z,z- 1}), 
r [0,i,i,o] = ([ w 5,W4,w 3 ],([z-l,y-l,x-l],{})), 
T [o,i,o,i] = (i w 5,yw 4 + W3y - W3,w 2 ],{[z,y - 1 + x],{2y - l,y - I}), 
T [o,i,o,o] = ([w5,w 4 ,w 2 ],{[z,y - l,x],{})), 
T [o,o,i] = ([-w 5 + w 4 ,W3 + 2zw 5 -w 5 ,w 2 -2zw 5 +w 5 ], 

([y,x- z],{z})), 
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Figure 3: Location of the sixth point for non null self-stress. 

T [o,o,o,i] = ([w 5 + 2yw4-W4,2w 3 y-W3 + w 5 ,W2 + w 5 ], 

([z,x-y],{2y-l})), 
T [o,o,o,o] = ([w5,w 3 -W4,w 2 ],{[z,2y - 1,2a;- 1],{ })), 

and the discriminant ideal is not principal: 

N = J = [y 2 z — yz, zx — z 2 , x 2 — y 2 — z 2 — x + y + z\. 

The generic solution is trivial (u/5 = W4 = W3 = W2 = 0). In this problem, the interesting 
non trivial solutions are given by the conditions over the parameters described by the variety 
of the discriminant ideal, which decomposes into 4 straight lines included in the hyperboloid 
x 2 — y 2 — z 2 — x + y + z = (illustrated in figure 3): 

Y(N) = Y(z, x-y)(J Y(y, x - z) (J Y(z, x + y - 1) (j V(y - 1, x - z). 

For this problem the Grobner basis wrt variables and parameters is already a comprehen- 
sive Grobner basis. 

6 Benchmarks 

For a set of examples taken from the literature we have applied the current implementation, 
release 2.3 in Maple 8 , of algorithm DISPGB using a 2 GHz Pentium 4 at 512 MB. Table 8 
summarizes the computing time of DISPGB, the total number of terminal vertices of the output 
tree, whether the discriminant ideal is principal or not, and whether the D-Grobner basis wrt 
>-xa is a CGB or not, joint by the number of failure cases for which it is not (0 if it is). The 
bases of the different examples are detailed below: 

• SI. [a(x + y), b(x + y),x 2 + ax}; 

• S2. [x\ , x\X2, x\x 2 , x\a + x 2 , X2X3 — x 2 , x 2 a, x|, x 2 a, a 2 ]; 

• S3, [x 3 — axy, x 2 y — 2y 2 + x]; 

• S4. [ax + y — 1, bx + y — 2, 2a; + ay, bx + ay + 1] ; 

• S5. [X4 — (04 — a 2 ), x 1 + X2 + X3 + X4 — (ai + 03 + 04), X1X3X4 — 040304, 

X1X3 + X1X4 + x 2 x 3 + X3X4 - (0404 + 0403 + 0304)]; 
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• S6. [vxy + ux 2 + x, uy 2 + x 2 ]; 

• S7. [y 2 — zxy + x 2 + z - 1, xy + z 2 — 1, y 2 + x 2 + z 2 — r 2 ]; 

• S8. [a — b + (xya — x 2 yb — 3a) 3 + (xyb — 3x6 — 56) 4 , xya — x 2 yb — 3a, 

xyb — 3x6 — 56]; 

• S9. [x + cy + bz + a, cx + y + az + 6, 6x + ay + z + c] ; 

• S10. See subsection 5.1; 

• Sll. [(d A d 3 R + r| - d 4 d 3 r 2 + d 2 d 2 - d A d\ - d\d 3 + d 4 d 3 + Z- R)t A 

+ {-2r 2 d A R + 2r 2 d\ + 1r 2 d A d\ - kr 2 d z d\ + 1r\ d A + 2r 2 d 4 )t 3 
- (2r 2 - 2ii + 4d 2 r 2 + 4d 2 + 2Z - 2d 2 A d 2 3 )t 2 
+ (-2r 2 d 4 R + 2r 2 d 4 dl + 2r 2 d 4 + 2r 2 d\ + 4r 2 d 3 d| + 2r\d A )t 
+ r\ + dfd 3 - d A d 3 R + d 4 d 3 r| + Z - R - d A d 3 + d\d\ + d 4 d|]; 

• S12. [a - / 3 c 3 - l 2 C!,b - l 3 s 3 - l 2 s u cf + s\ - 1, c§ + s§ - 1]; 

• S13. [ax 2 y + a + 36 2 , a(6 — c)xy + a6x + 5c]; 

• S14. [t 3 — cut 2 — -ui; 2 — utf 2 , t 3 — cvt 2 — vu 2 — vw 2 ,t 3 - cwt 2 — wu 2 — wv 2 }; 

• S15. [a + dsi, b — dc\,l 2 c 2 + l 3 c 3 - d, l 2 s 2 + / 3 s 3 - cs 2 + cf - 1, s| + c 2 - 1, 

s§ + c§-l]; 

• S16. See subsection 5.2. 

We have tested several other problems and in some of them only partial results have been 
reached. We detail two significative examples: 

• S17. [axt 2 + bytz - x(x 2 + cy 2 + dz 2 ), ayt 2 + bzxt - x(y 2 + cz 2 + dx 2 ), 

azt 2 + bxyt — x(z 2 + cx 2 + dy 2 )] 

• S18. [(3x 2 + 9v 2 -3v- 3x)tjtl + (3v 2 - 3v + 6vx - 3x + 3x 2 )^ 

+ (3v + 3v 2 + 3x 2 - 3x - 6vx)4 - 24v 2 ht 2 + 9v 2 - 3x + 3x 2 + 3u, 

(3x 2 + 9v 2 -3v- 3x)t 2 2 tl + (3v + 3v 2 + 3x 2 - 3x - 6vx)4 

+ (3f 2 - 3v + 6vx - 3x + 3x 2 )tl - 24v 2 t 2 t 3 + 9v 2 - 3x + 3x 2 + 3v, 

(3x 2 + 9v 2 -3v- 3x)tlt\ + (3f 2 - 3v + 6vx - 3x + 3x 2 )t\ 

+ {3v + 3v 2 + 3x 2 - 3x - 6vx)4 - 24v 2 t 3 h + 9v 2 - 3x + 3x 2 + 3v] 

For S17 [GoTrZaOO], DISPGB gets bogged down after computing 35 terminal vertices in 1375 
sec. It has been unable to finish the tree, and so neither rebuilding with the discriminantideal 
nor reducing the tree can have been achieved. The label of the 35th vertex is [1, 1, 0, 1, 0, 0], 
thus all vertices beginning with [0, 0, . . . have been already determined (the tree is built up 
in pre-order beginning with the vertices). 

S18 corresponds to the benzene molecule studied in [Em99]. The situation is similar to 
S17, getting bogged down after 45 seconds when the 9th vertex labelled [1, 1,0,0]] has been 
computed. 
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Identification 


CPU time 


# Final 


Discriminant 


Is CGB? 




(seconds) 


vertices 


is principal? 


(# failures) 


SI [We03] 


0.8 


2 


N 


Y (0) 


S2. [Gi87] 


1.2 


2 


Y 


N (1) 


S3. [Gom02, Du95] 


1.5 


2 


Y 


Y (0) 


S4. 


1.6 


2 


N 


Y (0) 


S5. [Kap95] 


1.6 


3 


Y 


N (1) 


S6. [Kap95] 


2.0 


4 


Y 


Y(0) 


S7. [Kap95] 


3.0 


2 


Y 


Y (0) 


S8. [SaSuNa03] 


4.4 


3 


Y 


Y (0) 


S9. Similar to [Si92] 


6.7 


10 


Y 


Y(0) 


S10. Subsection 5.1 


7.9 


7 


Y 


N (2) 


Simple robot 










Sll. [Co04] 


8.0 


6 


Y 


Y f0) 


Singular points 










S12. [RyOO] 

Rychlik robot 


8.2 


11 


Y 


N(l) 


S13. [SaSu03] 


8.2 


10 


Y 


Y(0) 


S14. [GoTrZaOO, De99] 
ROMIN robot 


9.6 


2 


Y 


N(l) 


S15. [GoRe93] 


18.2 


17 


Y 


N(2) 


S16. [GuOr04] 


21.3 


8 


N 


Y(0) 


Subsection 5.2 











Table 8: 
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